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ABSTRACT 

Studies of GRB host galaxies are crucial to understanding GRBs. However, 
since they are identified by the superposition in the plane of the sky of a GRB 
afterglow and a galaxy there is always a possibility that an association represents 
a chance alignment, rather than a physical connection. We examine a uniform 
sample of 72 GRB fields to explore the probability of chance superpositions. 
There is typically a ~ 1% chance that an optical afterglow will coincide with 
a galaxy by chance. While spurious host galaxy detections will, therefore, be 
rare, the possibility must be considered when examining individual GRB/host 
galaxy examples. It is also tempting to use the large and uniform collection of 
X-ray afterglow positions to search for GRB-associated galaxies. However, we 
find that approximately half of the 14 superpositions in our sample are likely to 
occur by chance, so in the case of GRBs localized only by an X-ray afterglow, 
even statistical studies are suspect. 

Subject headings: gamma rays: bursts 



Introduction 



From the earliest associations between gamma-ray bursts (GRBs) and host galaxies, 
burst en vironments have revealed many important clues to the nature of GRB progen i- 



tors (e.g. ISokolov et al.l l200ll : Bloom et all 12002 ; iFvnbo et al.ll2003l : iLe Floc'h et all 120031 ). 



As Swift accurately localizes hundreds of GRBs, statistical studies of GRB host galax- 
ies are becoming increasingly critical to under standing GRB formation mechanisms (e.g. 
Christensen et al.l 120041 ; IWainwright et al.l 120071 ). Such studies require that a large number 
of GRB host galaxies be identified. 
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A GRB host galaxy is identified by the superposition in the plane of the sky of the GRB 
afterglow (AG) and a galaxy. When an optical afterglow is detected, afterglow localization 
is precise to significantly less than an arcsecond. The superposition between a galaxy and 
an optical afterglow is, therefore, generally taken as proof that the galaxy in question is 
the host galaxy of the GRB. The fact that galaxies have no clear "edge" complicates this 
issue because there is no exact projected distance at which a GRB and a galaxy become 
associated. If only an X-ray afterglow is detected, the GRB can only be localized to within 
an error region of radius ~ 2 arcseconds. If a galaxy is detected within this error region it is 
generally assumed to be the host galaxy of the GRB, although confusion can arise when two 
or more sources are detected within or around a single X-ray afterglow error region. The 
issue of projected distance between a galaxy and the X-ray position also remains. 

With only visual associations between afterglows and galaxies used to identify GRB 
host galaxies, the possibility always exists for an incorrect association. Some afterglow and 
galaxy associations may be chance superpositions, with the galaxy either a foreground or 
background galaxy that is not physically associated with the GRB. These mistaken asso- 
ciations could cause confusion when analyzing data on GRB environments. When dealing 
with a statistical sample of GRBs, rather than a single case-study, results will only be im- 
pacted if a relatively significant number of the host galaxies have been mistakenly identified. 
Understanding the likelihood of galaxy mis- identification in a given sample of GRBs is, there- 
fore, imperative when conducting such studies. However, even a single mistaken association 
might generate confusion if the combination of GRB and galaxy characteristics are anoma- 



lous. For example, the potentially paradigm-shifting object GRB 06061 4 (e.g. ICobb et al. 



20061 : iDella Valle et al.ll2006l : iGal-Yam et al.ll2006l : ISchaefer & Xiaoll2006f ). a burst identified 



as low-redshift for which no SN was detected, could be a typical object of no special interest 
if the burst's purported host galaxy is actually just a random galaxy along the line-of-sight 
to the GRB. 

In this paper we investigate the fields of 72 GRBs in order to examine the general 
issue of associations between GRBs and host galaxies. In §2, we describe our data and 
analysis techniques. In §3, we address optical afterglow and host galaxy superpositions. In 
§4, localizations based on X-ray afterglows are considered. Afterglow-associated galaxies are 
compared to field galaxies in §5. We discuss our results in §6 along with some strategies for 
recognizing false hosts and conclude in §7. 
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Data 



Our optical data was obtained using the ANDICAM instrument mounted on the 1.3m 
telescope at Cerro Tololo Inter-American ObservatoryE] This telescope is operated as part 
of the Small and Moderate Aperture Research Telescope System (SMARTS) consortium^ 



2.1. GRB Selection 

The 72 bursts included in these analyses were, in general, selected based on their observ- 
ability from CTIO, limited to those bursts with declination < +35°. Observations were also 
limited by each burst's right ascension because all observations began within days post-burst. 
Bursts were occasionally not observed due to telescope scheduling time limitations. Both 
long-soft bursts (65) and short-hard bursts (7) were observed. Due to large and uncertain 
values of galactic reddening near the galactic plane, no GRBs with galactic latitude |6| < 10° 
are included in this sample. 

The redshifts, or redshift upper limits, of 27 bursts in our sample have been determined 
through either spectroscopic optical afterglow observations or observations of possible host 
galaxies. For this sample of bursts, redshifts range from 0.125 to 6.6. 

The GRBs included in this sample are listed in Table 1. 



2.2. Observations and Data Reduction 



The nightly data set for each burst consisted of 6 individual 360-second I-band obser- 
vations, taken at slightly offset telescope positions. Each individual image was reduced in 
the typical manner, with bias an d dark subtraction, flat fielding, and cosmic ray removal 
using the L.A. Cosmic program^ ( van Dokkum 2001 ). The individual images were median 
combined with maximum pixel rejection to produce a source-free initial fringe correction 
image. This background image was scaled and subtracted from the individual images, which 
were then aligned and combined. This combined image was used to produce a source mask 
and a new fringe correction image was produced by recombining the original images - only 



1 http:/ /www. astronomy.ohio-state.edu/ANDICAM 

2 http: / /www. astro. yale.edu/smarts 

3 htt p : / / www. astro .yale . edu/dokkum / lacosmic / 
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this time masking the sources in each image, rather than using minmax rejection. This new 
background was scaled and subtracted from the individual images, which were then aligned 
and combined to produce the final nightly science frame. 

For each GRB, between 1 and 11 usable nightly images were obtained (the mode being 
4 images). Particularly shallow images or nights with relatively poor seeing were excluded. 
Some images taken at early times post-burst were also excluded as they contained optical 
afterglow light. All usable nightly images were aligned and combined to produce a final deep 
frame for each GRB. The frame edges were cropped so that each final frame is of uniform 
depth over the entire field. Field size varied slightly but is typically 5'6 x 5'A. A few fields 
contained saturated stars with diffraction spikes; the area immediately surrounding these 
stars is excluding from all of the following analyses. 

Secondary standard stars in each field were ph otometrically calibrated by comparison, 
on photometric nights, with Landolt standard stars (lLandoltlll992l ). Photometric calibration 
is typically accurate to 0.05 magnitudes. No photometric observations were obtained for 
4 GRBs and those fields are calibrated using USNO-B1.0 1 2 magnitudes, with typical cali- 
bration errors of 0.2 magnitudes. The reddening corrected (ISchlegel et al.lll998l ) 3a limiting 
magnitudes for point sources in these images ranges from 21.2 mag to 23.3 mag, with a 
median of 22.4 mag. 

The fields are astrometrically calibrated using USNO-B1.0 stars, with statistical error 
of < 0.2". The CCD pixel scale is 0.37" /pixel, so that the astrometry is accurate to within 
a single pix el. X-r a y afte rglow coordinates and 90% confidence error radii for all bursts are 
taken from iButlerl (120071 ). Error region radii vary from / .'5 to 6'.'9, with a median of l'!9. 



The data are summarized in Table 1. 



2.3. Galaxy Detection 



All objects in the field were cataloged using SExtractor (IBertin fc Arnoutsl 119961 ) with 
a 2a detection threshold. Using SExtractor's neural network star/galaxy classifier, objects 
with CLASS_STAR> 0.8 are considered stars and not included in the following analyses. In 
crowded starfields, SExtractor's ability to differentiate between stars and galaxies is reduced. 
This effect is minimized in our analyses, however, as none of the fields considered here are 
within 10° of the galactic plane. Galaxy I-band magnitudes are given by MAG_AUTO, with 
zeropoints appropriately adjusted for each field. These magnitudes are generally accurate to 
within a few tenths of a magnitude. 
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Only galaxies detected by SExtractor are considered in these analyses. At the position 
of a few afterglows there are low significance galaxies that are not identified by SExtractor. 
For consistency, these visual identifications are not included. 

Since our observations are limited to galaxies brighter than ~ 23 mag, we would not ex- 
pect to be able to detect all the potential host galaxies in our images. Galaxies go undetected 
primarily because of their large redshifts, though dwarf galaxi es might be missed ey en at low 
redshifts. Examination of the VVDS-CDFS galaxy catalog (ILe Fevre et al.ll2004l ) suggests 
that nearly all of the detected galaxies in these fields are likely to be at a redshift lower than 
z = 1.5 and about 75% have redshift z < 0.7. 

To understand our ability to detect GRB host galaxies, it would be useful to know the 
characteristics of a typical GRB host. In pie-Swift GRB samples, however, galaxies associ- 



ated with GRBs have been noted to be a rather heterogeneous group (e.g. IConselice et al. 



20051 : ILe Floc'h et al. 



2003 ) covering a wide range of galaxy types and absolute magni tudes 



16 to —21) ( Sokolov et al. 



Conselice et al. \ 


>005; 


Wainwright et al. 


2007 


). As- 


2 ( 


Fukugita et al 


199, 


5), the brightest host galaxies 



suming a typical galaxy color of B — 1 
could have have / = —23 and might be detected in our images to z ~ 1. Observed galaxy 
brightness would, of course, depend on the necessary fc-correction, which depends strongly 
on galaxy morphology. The star-forming galaxies of long-duration bursts would be favored 
for detection at higher redshifts because of their strong rest-frame UV emission. Unfortu- 
nately, redshift limits for detecting the dimmest host galaxies are much more severe. These 
dwarf hosts may be common in the local universe, but many could not be detected in our 
images to even moderate redshifts given our magnitude limits. 



3. Optical Afterglow and Galaxy Coincidence 

The position of optical afterglows can be determined to within sub-pixel accuracy and, 
therefore, can significantly constrain the probability of a GRB/galaxy association. The 
probability of the optical afterglow falling on any galaxy-covered pixels is determined by 
dividing the number of pixels covered by galaxies in each field by the total number of pixels 
in that field. 



3.1. ISOAREAFJMAGE area 



One measure of the pixel area of each galaxy is given by the isophotal SExtractor 
parameter ISOAREAFJMAGE. This value disregards the lower significance outskirts of 
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each galaxy, so it is a conservative estimate of galaxy area. If an optical afterglow were to 
fall within the ISOAREAFJMAGE area of a galaxy, that galaxy would be regarded as a 
strong host-galaxy candidate. This excludes cases in which GRBs fall at the outskirts of 
galaxies - which are bound to occur. Hence these figures serve strictly as a lower limit. 

The ISOAREAFJMAGE area covered by galaxies in each field ranges from 0.2% to 
3.5%, with a median of 1.0% (see Figure 1 and Table 1). Summing the overlap probability 
over all 72 fields yields an expectation of 0.86 observed galaxy/optical afterglow coincidences 
- including even those fields for which no optical afterglow detection was made. For the 47 
fields in which optical afterglow was detected from either short- or long-duration bursts, 
there is an expectation of 0.54 observed galaxy/optical afterglow coincidences. In fact, this 
sample contains 7 such coincidences (GRBs 050416a, 050724, 050826, 060505, 060614, 061021 
and 061121). The chance of having randomly observed that many coincidences is less than 
0.004%. 

Since the total galaxy-covered area in each field depends on limiting magnitude, which 
varies from field to field, it is useful to consider a more homogeneous selection of galaxies. 
As deeper and deeper observations are obtained, more and more galaxies will be detected. 
Therefore, deeper imaging will serve to increase the chance for spurious associations. By 
imposing a magnitude cutoff, therefore, we produce a minimum value for the possibility of 
chance superposition. For the homogeneous sample, we limit our analysis to galaxies with 
Galactic extinction corrected magnitudes brighter than / = 21.5. This value is determined 
by producing a histogram of all the galaxies being used in this analysis, with bin size of half a 
magnitude, and selecting the midvalue of the bin containing the maximum number of galaxies 
(see Figure 2). Beyond this magnitude, our galaxy sample begins to become significantly 
affected by incompleteness. In individual fields, this galaxy completeness "turnover" value 
ranges from 20 to 22, with a median and mode of 21.5. In the fields that are incomplete 
at 21.5 mag, there will be "missing" galaxies. This means that the probability of chance 
superposition with a / < 21.5 galaxy will be slightly underestimated. 

Not surprisingly, the exclusion of all / > 21.5 galaxies reduces the field coverage by a 
small amount, so that the median galaxy-covered area is only reduced to 0.9%. Of the 7 
GRB/galaxy coincidences noted above, 4 are with I < 21.5 galaxies (GRBs 050724, 050826, 
060505, 061021) The expected number of observed coincidences with / < 21.5 galaxies is 0.79 
(or 0.50 if considering only the 47 GRBs with detected optical afterglows), so the chance 
of having observed 4 at random is < 1%. Even in this reduced sample, there remains a 
meaningful overdensity of galaxy/optical afterglow coincidences. However, it is plausible 
that one or more of these associations could be a coincidence. 
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3.2. Ellipse Area 

Without fully understanding how all GRBs are produced, the exact placement of GRBs 
within their galaxies cannot be accurately predicted. Long-duration GRBs, for example, 
may require the kind of rapid star formation regions often found in galactic spiral arms, 
while short-duration GRBs may favor the fringes of galaxies if they are formed by the 
mergers of compact remnants. We have, therefore, been somewhat too restrictive in as- 
suming that a galaxy will only be identified as a host if the GRB occurs within the galaxy's 
ISOAREAFJMAGE area. 

To consider the more general situation of a optical afterglow at a given observed position 
relative to a nearby galaxy, we define galaxy area as an ellipse having major and minor axes 
of length ux AJMAGE and ux BJMAGE, where AJMAGE and BJMAGE are SExtractor 
parameters and v is a simple scaling factor. In general, an ellipse with v — 3 is visually 
coincident with the extent of the galaxy it describes. When the shape of each galaxy is 
defined in this way, galaxy area is equal to v 2 nx AJMAGE x BJMAGE. 

In Figure 3, we plot the probability of chance alignment for an afterglow contained within 
a galaxy ellipse with scale factor v. The probability of chance alignment rises quadratically 
with v until v is so large that the galaxies significantly overlap with one another. We plot 
both the entire sample of galaxies and only those galaxies with I < 21.5. Excluding the 
dimmer galaxies only changes the random overlap probability by a few tenths of a percent 
at large v. On the right hand side of the graph, the number of expected chance coincidences 
in our sample of 72 GRBs is shown. 

If the scale factor is set to v = 3 so that galaxies are described by ellipses which are 
visually coincident with the extent of each galaxy, then the ellipse area covered by galaxies 
in each field ranges from 0.4% to 5.2%, with a median of 1.6%. Summing the overlap 
probability over all 72 fields yields an expectation of 1.3 observed galaxy/optical afterglow 
coincidences. For the 47 fields in which optical afterglow was detected there is an expectation 
of 0.82 observed galaxy /optical afterglow coincidences. The chance of having observed the 
7 coincidences seen in this sample is less than 0.04%. Limiting ourselves to galaxies with 
/ < 21.5, the median galaxy-covered is 1.4%. The expected number of observed coincidences 
with / < 21.5 galaxies is 1.2 (or 0.73 if considering only the 47 GRBs with detected optical 
afterglows) so the chance of having observed 4 at random is < 3%. 
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3.3. Comparing Short- and Long-Duration GRB Fields 



A comparison of short- and long-duration fields is of interest, although the sample of 
short duration bursts in this sample is limited to only 7 fields (versus 65 long-duration 
bursts). The galaxy area coverage spans a similar range of values between the two samples. 
The median area (ISOAREAFJMAGE area) covered by all galaxies in short-GRB fields 
(~ 0.009) is only slightly lower than the median area covered in long-GRB fields (~ 0.01). 
Since short bursts may occur further outside their galaxies than long-duration bursts, the 
galaxy area over which a true association would be physically plausible may be larger for 
short bursts than for long bursts. 

In this sample, three short-GRBs had detected optical afterglows, but there is only 
one short-burst galaxy /optical afterglow association. This association frequency (1 of 3) is 
comparable to that of the long-GRBs (6 of 44). If short bursts are assoc iated with early-type 



galaxies, they may be preferentially located in local galaxy clusters (e.g. lPedersen et al.ll2005 



Bloom et al.ll2006l ). If so, a significantly higher frequency of associations might be expected 
for short bursts than long bursts. As this is not yet confirmed, these observations are still 
consiste nt with recent eviden ce that suggests that not all short bursts are limited to local 
clusters (IBerger et al.ll2007aU bl). 



4. X-ray Afterglow and Galaxy Coincidence 

For a non- negligible population of Swift GRBs (nearly 50%), no optical afterglow is 
detected. Underluminous optical afterglows and heavy line-of-sight extinction may account 
for many "dark" bursts. An optical afterglow might also go undetected due to observing 
constraints, such as the relative position of the sun or the moon to the burst's coordinates or 
poor weather conditions at optimum observing sites. Regardless of the reason, such a non- 
detection generally means that the burst will only by localized to within a few arcseconds by 
XRT observations of the burst's X-ray afterglow. XRT observations comprise a dataset that 
is significantly more uniform than deep, ground-based optical afterglow observations that, by 
necessity, are obtained from a large number of different instruments at varying times post- 
burst. This homogeneity, combined with the fact that the X-ray localized dataset contains 
nearly double the number of optically localized GRBs, makes it tempting to analyze GRB 
host galaxies based exclusively on X-ray AG positions. We, therefore, examine our data to 
determine how significantly such an analysis might be impacted by the presence of falsely 
identified hosts. 

Each X-ray afterglow error region of the GRBs in this sample was examined for coin- 
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cident galaxies, and 14 are found to overlap with one or more SExtracted galaxiesjj Any 
overlap between the circle defined by the X-ray afterglow error region and a galaxy ellipse 
with v = 3 is considered to be a coincidence. This definition of galaxy area is used instead of 
the ISOAREAFJMAGE area because an ellipse defines a clear "edge" to the galaxy while 
SExtractor does not output the exact boundaries of a galaxy. 

We then investigated the probability that one or more galaxies will fall within any 
random region that is the same size as a burst's X-ray afterglow error region. For each 
GRB field, 500 random positions were selected and each position was assigned a region 
with a radius equivalent to the X-ray afterglow error radius of the corresponding GRB. All 
region/galaxy overlaps were then counted. Only a handful of overlaps occur in some fields, 
while in others more than a quarter of all random regions contain a galaxy. The median 
overlap probability is 6.9%. The probability of one or more galaxies overlapping any region 
the size of the burst's X-ray afterglow error region is given in column 8 of Table 1. The 
overlap probability in each field strongly depends on both the radius of the random regions 
and the galaxy population density. The burst with the largest X-ray error region is GRB 
050412, with a radius of 6'.'9. In this field, any randomly placed region of that size has a 
32% chance of overlapping with a galaxy. Note that the density of the galaxy population 
is dependent on image depth, as shallower images will contain fewer dim galaxies. Chance 
probabilities will, therefore, increase with image depth. Figure 4 shows a histogram of the 
probability of region/galaxy overlap in the 72 GRB fields examined. Arrows on the graph 
indicate the fields in which the X-ray afterglow error region did coincide with one or more 
galaxies. 

The expected number of observed overlaps created by chance in this GRB sample is 6.7, 
which is obtained by summing the overlap probability over all fields. This is significantly less 
than the 14 actually observed. Using a Monte Carlo simulation (see Figure 5), the probability 
of having observed 14 overlaps in this sample is found to be only 0.5%. Thus, there is a clear 
over-density of galaxies in GRB X-ray afterglow error regions. But the galaxies associated 
with the X-ray afterglow error regions do not have to be the true host galaxies of the GRBs. 
In fact, of the X-ray afterglow error regions in this sample that contain a galaxy, it is likely 
that approximately half contain a galaxy that is not associated with the GRB. 

Approximately 250 Swift bursts have been detected and well-localized by their X-ray 
afterglows using the XRT. Based simply on the median probability of overlap in these fields, 
there is only a 10 -8 chance that no "false hosts" have been detected. In fact, if our sample 
represents a typical distribution then there could be over 20 such detections in the entire Swift 



4 The X-ray afterglow error region of GRB 060505 overlaps with 2 galaxies. 
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sample. Improving the limiting magnitudes in each field would only serve to significantly 
increase the number of galaxy detections. Indeed, as image depth increases the detection 
of one or several galaxies within a given X-ray afterglow error region is inevitable. Clearly, 
caution is required when identifying a galaxy in a X-ray afterglow error region as the host 
galaxy of the GRB. Note, in particular, that relatively few optical afterglow have been asso- 
ciated with short-duration bursts, yet many claims as to the nature of short- duration bursts 
have recently been made on the basis o f their X-ray afterglow and galaxy associations (e.g. 



Pedersen et al.ll2005l ; iBloom et al.ll2006l ). The likelihood of misidentified hosts is, therefore, 



of particular concern when considering short bursts. 



5. Field Galaxies versus Afterglow- Associated Galaxies 

Having collected the observable properties of all the galaxies in these 72 fields, we can 
determine if the galaxies associated with either X-ray or optical afterglows represent an 
anomalous sample of the field galaxies. The observable properties considered are magnitude, 
pixel area (either ISOAREAF_IMAGE area or ellipse area) and ellipticity (e), where e = 1- 
(BJMAGE/AJMAGE). These observable properties are, of course, a function of intrinsic 
size, shape, brightness and redshift. 

Figure 6a shows galaxy ellipse area, with v — 3, versus magnitude. Monte Carlo sim- 
ulations show that the galaxies associated with X-ray afterglow error regions are a typical 
sample of all the galaxies in the field over this dimension. This result is not surprising, since 
nearly half these galaxies could actually be just field galaxies, rather than GRB-associated 
galaxies. The galaxies associated with optical afterglow, however, do appear slightly different 
than the field galaxies at the ~ 2a level. This seems to be due to the fact that these galax- 
ies are somewhat brighter than typical field galaxies: ~ 30% of optical-afterglow-associated 
galaxies have I < 18 mag, while this is true of only ~ 3% of field galaxies. To test this, we 
split the entire sample of galaxies into three magnitude bins: bright: / < 18, intermediate: 
18 < / < 22 and faint: I > 22. Of the 7 galaxies associated with optical afterglow, 2 are 
bright, 3 are intermediate and 2 are faint. We then repeat the Monte Carlo simulation but 
require that for each run the 7 galaxies randomly picked from the entire sample include 2 
bright, 3 intermediate and 2 faint galaxies. With the addition of this magnitude selection, the 
optical afterglow-associated galaxy sample becomes indistinguishable from the field galaxies. 

Similar analyses are performed for galaxy ISOAREAFJMAGE area versus magnitude, 
ellipticity versus magnitude, galaxy ellipse area versus ellipticity and galaxy ISOAREAFJMAGE 
area versus ellipticity (see Figure 6b-e). The X-ray afterglow error regions are a typical sam- 
ple of all the galaxies in the field over these dimensions. The optical afterglow-associated 
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galaxies are atypical on the ~ 2a level when comparing ISOAREAFJMAGE area versus 
magnitude and ellipticity versus magnitude. Again, this seems to be due to the fact that 
these galaxies are somewhat brighter than the average field galaxies. Since 6 out of 7 of 
these galaxies are associated with long-duration GRBs, the fact that long-duration GRB 
host generally have strong ongoing star-formation may account for their relative brightness 
when compared with the field galaxies in general. 

6. Discussion 

Unlike the fairly uniform images considered here, most observations of GRB fields are 
obtained with a wide range of depth and resolution. Improving depth and resolution will 
result in the detection of more galaxies, in which case the values calculated here serve only 
as a lower limit to the probability of chance association. While increased resolution may 
seem to make the "edge" of a galaxy more defined, we would caution that knowing the exact 
optical extent of the galaxy is not helpful for identifying a GRB's host without first fully 
understanding how GRBs trace light. 

An additional complication is lensing, which will always occur when a background GRB 
is observed through a foreground galaxy. Lensing will increase the brightness of the GRB and 
its afterglow and shift the observed position of the GRB's afterglow relative to the lensing 
galaxy. In principle, this could increase the probability of observing chance superpositions. 
However, we demonstrate in the Appendix that this effect will be negligible in most cases. 

Currently, there are only a few situations in which any kind of mis-match between a 
GRB and proposed host is likely to be noted. One example of this would be when no SN is 
detected in a low-redshift galaxy associated with a long-duration GRB. A second example 
would be the association of a long-duration burst with an elliptical galaxy with little to no 
ongoing star- formation. In either event, it is difficult to make a firm conclusion because of 
the possibility that not all GRBs are produced in the canonical fashion. One can question 
the cause of a GRB as easily as question its redshift. 

To avoid contaminating host galaxy samples with false-hosts, it would be useful to have 
a way to separate true and false hosts. Unfortunately, no exact method of doing so exists. 
Some strategies are: 

• Gamma-ray luminosity indicators: If luminosity indicators are both correct and 
universal, they can be used to calculate a burst's redshift. However, a discrepancy could 
either indicate a chance superposition or a very unusual GRB, so that interpretation 
of such a mismatch is unclear. Lensing effects could also invalidate this method. 
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Detection (or non-detection) of SNe: This is only meaningful if long-duration 
bursts are all formed in core-collapse SNe and are similar in brightness to GRB 
980425/SN 1998bw, which is assumed to be the archetypal GRB-related SN. This 
is further limited to GRBs that occur at relatively low redshifts (z < 0.7). Host-galaxy 
extinction is also a complicating factor that may obscure an underlying SN. 

Classification of associated galaxy: Long-duration GRBs should be found in star- 
forming, late-type galaxies, while short-duration GRBs should be found in quiescent 
early-type galaxies. This assumes, however, that all GRBs are produced in the canon- 
ical fashion. 

Optical afterglow spectral absorption features: A GRB must occur behind 
the absorbing material, so line detection actually only gives a lower limit on redshift 
(though this is generally taken to be the GRB's redshift). Redshifts derived in this 
manner can be compared to redshifts measurements of the proposed host galaxy after 
the afterglow has faded. Only a very bright afterglow could possibly show absorption 
features independent of the superimposed galaxy, however. Alternatively, an afterglow 
with no absorption features might be assumed to be at low redshift because of the 
complete lack of absorption features along the line-of-sight. 

Visual detection of true host galaxy: In deep images, the true host galaxy may 
be detectable through the foreground galaxy. An example of this situation is GRB 
060912a which occurred in a z = 0.937 ga laxy that overlaps with the outskirts of a 



large z = 0.0936 galaxy (ILevan et al.l 120071 ). 



X-ray afterglow spectral emission lines: Such emission lines have been suggested 
to be produced by material associated with the GRB (outflows, disks, etc.) so they 
could serve as accurate redshift indicators. This is very unproven method, however, 
because line identification is inexact, reported lines h ave been at only low significance, 



and few GRB X-ray spectra contain these lines (e.g. ISako et al.l 120051 1 . 



Strong lensing of GRB afterglows: This would be extremely unlikely for a single 
foreground galaxy, (as shown in Appendix A) but could potentially be caused by a 
foreground galaxy cluster. 



7. Conclusions 



We have studied the galaxy population surrounding a sample of 72 GRBs. Typically 
1% of the sky near the positions of GRBs is covered by galaxies with / < 21.5. With ~ 125 
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Swift GRBs with detected optical afterglows, the probability that no chance alignments have 
been detected is < 28%. Indeed, approximately 1 superposition between the GRB and a 
foreground galaxy is expected. While it is possible that no such chance alignments have yet 
been observed, as Swift detects more and more GRBs with optical afterglows the likelihood 
of such an event only increases. While most GRB/galaxy associations noted in the literature 
are almost certainly correct, caution is required when making sweeping conclusions from 
only one or two GRB/host galaxy associations. 

We have also considered galaxies associated with X-ray afterglows. Over 250 X-ray 
afterglows, with typical error radii of ~ 2", have been detected by Swift. These numbers 
guarantee that some galaxies associated with X-ray afterglow will be falsely identified as 
hosts. In fact, approximately half of the 14 X-ray afterglow error region/galaxy coincidences 
in our sample may exist only by chance. Even with large samples, therefore, using X-ray 
afterglow-identified host galaxies to draw conclusions about GRBs may cause confusion. 

We thank SMARTS observers D. Gonzalez, J. Espinoza and A. Pasten for their dedica- 
tion and S. Tourtellotte for assistance with optical data reduction. We are also grateful to 
P. G. van Dokkum and P. Natarajan for useful discussions about galaxies and lensing. This 
work is supported by NSF Graduate Fellowship DGE0202738 to BEC and NSF/AST grant 
0407063 and Swift grant NNG05GM63G to CDB. 



A. Lensing 

Though the chance of a background GRB coincidently falling behind an unassociated 
galaxy is only about 1%, if a chance alignment does occur then the GRB will be lensed. The 
magnification produced by this lensing, of course, depends strongly on galaxy mass, observer 
to lens distance, lens to source distance, and the extent of the alignment. This lensing effect 
actually increases the likelihood of detecting a chance association because it increases the 
observed brightness of a GRB's optical afterglow. Host galaxy/optical afterglow associations, 
by definition, require the detection of an optical afterglow, and brighter optical afterglows 
are more likely to be detected than dim afterglows. 

Lensing presents an additional complication because lensed GRBs might not follow 



expected luminosity relationships. The Amati relationship (lAmati et al.ll2002l ; lAmatill2006l ). 
for example, correlates gamma-ray spectral peak energy with a burst's isotropic equivalent 
radiated energy. From measurements of peak spectral energy and gamma-ray flux, therefore, 
a GRB's redshift can be estimated. Lensing should not alter the gamma-ray burst spectral 
shape so that peak energy remains unchanged. The increase in gamma-ray flux produced by 
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lensing, however, would result in an underestimated burst distance (z oc luminosity /flux). 
This is problematic for detecting inconsistent GRB and galaxy redshifts because lensing will 
move the pseudo-redshift of the GRB toward the redshift of the foreground galaxy. 

Treating the lensing galaxy as a single isothermal sphere, the strongest lensing effect 
would be produced by a low-redshift galaxy lying directly along the line-of-sight to a high- 
redshift GRB (angular offset between source and lens < 0'.'3). Random chance, however, 
favors less exact alignments. At greater angular separations (> 0'.'3, where magnification 
only depends weakly on angular separation), only very large galaxies (>L* or a > lOOkm/s) 
can produce more than a factor of 2 in magnification of the source. Consider a scenario in 
which a 0.1L* galaxy at z=0.3 lenses a GRB at z=3, with an angular separation of Of! 5. This 
produces only a 15% increase in brightness. Lensing, therefore, only becomes a concern when 
a burst aligns very closely with a >L* galaxy. However, there is only a 10~ 4 chance that the 
line-of-sight to a GRB at redshift z = 5 would pass within 0f!3 of the gravitational center of a 
>L* galaxy. For bursts that occur at lower redshift, this probability only decreases. Hence, 
while lensing will always occur with chance alignments, the effects will generally be of little 
consequence. 
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Table I. Properties of 72 SMARTS /-band GRB fields 



X-ray AG Coordinates and Radius a Optical Limiting AG Overlap Probability Observed Galaxy /AG Overlaps 



\jSXD 






I ) 




/ magk 


cpticai 


X-ray 


(and / mag of associated galaxy) 


050128 


14:38:17.66 


-34:45:54.7 


2.7 




22.4 


0.012 


0.118 




050219a 


11:05:39.13 


-40:41:02.6 


3.5 




22.5 


0.010 


0.154 


X-ray (18.9) 


050223 


18:05:32.66 


-62:28:20.4 


3.2 




22.3 


0.017 


0.194 


X-ray (21.1) 


050315 


20:25:54.21 


-42:36:02.5 


1.4 


yes 


22.5 


0.018 


0.078 




050318 


03:18:51.04 


-46:23:43.5 


2.7 


yes 


22.4 


0.010 


0.088 




050401 


16:31:28.84 


+02:11:14.5 


1.8 


yes 


22.6 


0.009 


0.088 




050408 


12:02:17.36 


+10:51:10.5 


1.2 


yes 


22.4 


0.011 


0.048 




050412 


12:04:25.18 


-01:12:00.8 


6.9 




22.6 


0.011 


0.322 


X-ray (21.7) 


050416a 


12:33:54.57 


+21:03:26.9 


0.6 


yes 


22.7 


0.012 


0.036 


X-ray + optical (22.2) 


050502b 


09:30:10.06 


+16:59:46.5 


1.0 


yes 


22.3 


0.006 


0.014 




050509b d 


12:36:14.06 


+28:59:07.2 


3.4 




22.3 


0.023 


0.190 




050520 


12:50:05.94 


+30:27:03.2 


2.4 




22.1 


0.009 


0.038 




050525 


18:32:32.67 


+26:20:21.2 


2.1 


yes 


21.7 


0.009 


0.044 




050714b 


11:18:47.75 


-15:32:49.3 


2.1 




22.0 


0.008 


0.062 


X-ray (20.9) 


050724 d 


16:24:44.32 


-27:32:26.4 


5.0 


yes 


21.6 


0.035 


0.572 


X-ray + optical (17.4) 


050726 


13:20:12.16 


-32:03:51.0 


3.7 




22.8° 


0.025 


0.208 




050730 


14:08:17.12 


-03:46:16.3 


1.4 


yes 


22.0 


0.005 


0.024 




050801 


13:36:35.51 


-21:55:42.7 


5.0 


yes 


22.1° 


0.011 


0.196 




050803 


23:22:37.87 


+05:47:09.8 


1.8 




22.5 


0.019 


0.090 




050822 


03:24:27.22 


-46:02:00.0 


0.7 




22.9 


0.018 


0.034 




050826 


05:51:01.69 


-02:38:37.6 


2.6 


yes 


21.5 


0.005 


0.056 


X-ray + optical (19.5) 


050827 


04:17:09.58 


+18:12:00.4 


1.3 




21.3 


0.005 


0.018 




050915a 


05:26:44.86 


-28:00:59.9 


1.4 


yes 


23.0 


0.015 


0.060 




050922b 


00:23:13.19 


-05:36:15.6 


2.7 




22.5 


0.008 


0.074 




051001 


23:23:48.73 


-31:31:23.3 


1.5 




23.1 


0.014 


0.050 




051006 


07:23:14.03 


+09:30:21.9 


4.3 


yes 


22.5 


0.028 


0.292 




051021a 


01:56:36.41 


+09:04:04.3 


2.1 


yes 


22.8 


0.007 


0.072 




051117b 


05:40:43.21 


-19:16:27.2 


2.0 




22.8 


0.010 


0.078 


X-ray (20.3) 


051210 d 


22:00:41.26 


-57:36:46.5 


2.9 




21.9 


0.009 


0.108 




060108 


09:48:02.00 


+31:55:08.0 


1.3 


yes 


22.4 


0.010 


0.046 




060115 


03:36:08.18 


+17:20:44.3 


2.6 


yes 


22.1 


0.005 


0.046 




060116 


05:38:46.14 


-05:26:15.2 


3.1 


yes 


22.3 


0.002 


0.028 




060204b 


14:07:15.05 


+27:40:36.9 


1.6 


yes 


22.4 


0.024 


0.068 




060211b 


05:00:17.12 


+14:56:58.1 


1.9 




21.3 


0.005 


0.040 




060213 


09:26:23.05 


-09:06:43.4 


3.2 




23.1 


0.021 


0.186 




060313 d 


04:26:28.41 


-10:50:40.7 


2.4 


yes 


22.6 


0.008 


0.066 




060505 


22:07:03.39 


-27:48:53.0 


2.5 


yes 


22.7 


0.010 


0.104 


X-ray + optical (17.5) and X-ray (20. 


060526 


15:31:18.23 


+00:17:05.4 


1.6 


yes 


22.1 


0.011 


0.044 




060602a 


09:58:16.93 


+00:18:12.8 


2.3 


yes 


21.9 


0.004 


0.018 




060604 


22:28:55.01 


-10:54:55.9 


1.0 


yes 


21. 2 C 


0.003 


0.008 




060614 


21:23:32.04 


-53:01:37.0 


1.4 


yes 


23.3 


0.022 


0.076 


X-ray + optical (21.9) 


060708 


00:31:13.79 


-33:45:33.3 


1.4 


yes 


22.5 


0.011 


0.056 




060714 


15:11:26.43 


-06:33:59.5 


1.8 


yes 


22.2 


0.009 


0.062 




060719 


01:13:43.69 


-48:22:51.0 


1.5 


yes 


22.9 


0.010 


0.066 




060729 


06:21:31.90 


-62:22:12.6 


0.5 


yes 


22.3 


0.009 


0.020 
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Table 1 — Continued 





X-ray AG Coordinates and Radius a 


Optical 


Limiting AG Overlap Probability 


Observed Galaxy/ AG Overlaps 


GRB 


RA DEC (") 


AG? 


/ mag b Optical c X-ray 


(and / mag of associated galaxy) 



060801 d 


14:12:01.35 


+16:58:53.7 


2. i 




22.1 


0.005 


0.036 


060805 


14:43:43.46 


+12:35:11.8 


1.8 




22.1 


0.009 


0.050 


060814 


14:45:21.29 


+20:35:10.7 


0.9 


ves 


21.6 


0.004 


0.008 


060904b 


03:52:50.45 


-00:43:30.0 


1.7 


yes 


22.2 


0.018 


0.080 


060912a 


00:21:08.11 


+20:58:18.9 


1.9 


yes 


22.5 


0.017 


0.088 


060919 


18:27:41.89 


-51:00:51.1 


2.6 




22.0 


0.022 


0.222 


060923c 


23:04:28.36 


+03:55:28.4 


1.8 


yes 


21.8 C 


0.013 


0.078 


060926 


17:35:43.86 


+13:02:19.0 


4.6 


yes 


21.2 


0.016 


0.250 


060927 


21:58:12.23 


+05:21:52.2 


3.8 


yes 


22.8 


0.012 


0.170 


061004 


06:31:10.93 


-45:54:23.8 


1.6 




22.2 


0.011 


0.068 


061006 


07:24:07.56 


-79:11:56.6 


2.4 


yes 


22.0 


0.005 


0.056 


061007 


03:05:19.55 


-50:30:01.8 


1.5 


yes 


22.4 


0.012 


0.042 


061021 


09:40:36.04 


-21:57:04.9 


1.3 


yes 


22.4 


0.010 


0.028 


061025 


20:03:37.18 


-48:14:28.9 


3.4 


yes 


22.0 


0.011 


0.144 


061102 


09:53:37.84 


-17:01:26.0 


2.9 




22.4 


0.010 


0.100 


061110a 


22:25:09.92 


-02:15:31.6 


1.6 


yes 


22.4 


0.020 


0.084 


061121 


09:48:54.50 


-13:11:43.0 


1.0 


yes 


22.6 


0.008 


0.026 


061201 d 


22:08:32.09 


-74:34:48.2 


3.6 


yes 


22.1 


0.006 


0.104 


061202 


07:02:05.54 


-74:41:54.6 


1.5 




22.2 


0.008 


0.052 


061217 d 


10:41:39.32 


-21:07:22.1 


3.8 




22.5 


0.009 


0.142 


070224 


11:56:06.57 


-13:19:49.4 


1.7 


yes 


22.7 


0.016 


0.076 


070306 


09:52:23.32 


+10:28:55.4 


1.1 


yes 


22.4 


0.008 


0.032 


070318 


03:13:56.79 


-42:56:47.7 


1.4 


yes 


21.8 


0.003 


0.002 


070330 


17:58:10.53 


-63:47:35.0 


2.4 


yes 


22.2 


0.018 


0.136 


070419b 


21:02:49.91 


-31:15:47.2 


1.7 


yes 


22.5 


0.021 


0.096 


070429a 


19:50:48.71 


-32:24:15.1 


1.9 




22.2 


0.020 


0.126 


070508 


20:51:12.10 


-78:23:03.4 


2.1 


yes 


22.2 


0.008 


0.070 



X- 



X- 



X 



X- 



ray (21.1) 
ray (22.0) 

ray + optical (21.4) 
■ray + optical (22.2) 



a X-ray afterglow coordinates and 90% confidence error radii for all bursts are taken from lButlerl j2007t) ■ 

b 3tr limiting / magnitude for point sources in the field, corrected for Galactic reddening. Unless otherwise indicated, photometric 
calibration is typically accurate to 0.05 magnitudes. 

c Based on the ISOAREAF JMAGE area of each galaxy. 

d Short-duration burst. 

°Photometrically calibrated using USNO-B1.0 12 magnitudes, typically accurate to 0.2 magnitudes. 
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Fig. 1. — Histogram of the probability of any pixel in a GRB field falling inside of a galaxy 
[solid line - median of 1.0%) or inside a galaxy with magnitude I < 21.5 (dotted line - median 
of 0.9%). The area of each galaxy is given by the SExtractor output ISOAREAFJMAGE. 
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Fig. 2. — Histogram of all galaxies in this sample. The magnitude cut of 1=21.5 is chosen as 
the midvalue of the bin containing the maximum number of galaxies. The sample of dimmer 
galaxies is significantly affected by incompleteness. 
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Fig. 3. — The probability of chance overlaps for (solid squares) all galaxies or (open circles) 
galaxies with I < 21.5 for increasing values of v, where v is the scale factor used in defining 
galaxy size. On the right is shown the corresponding expectation value for the 72 GRBs 
in this sample. The v values of the 7 GRBs with optical afterglow/galaxy coincidences are 
marked. 
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Fig. 4. — Histogram of the probability of a random X-ray error region/galaxy overlap in the 
72 GRB fields examined. The median probability is 6.9%. The arrows indicate the actual 
probability of each field in which the X-ray afterglow error region did coincide with one 
or more galaxies. The arrows over > 0.3 indicate fields having actual probabilities of 32% 
and 57%. A Kolmogorov-Smirnov (K-S) test on these data gives P ~ 0.4, which indicates 
that the distribution of probabilities in the 14 fields containing X-ray afterglow-associated 
galaxies is not significantly different from the distribution of the entire sample. 
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Fig. 5. — Monte Carlo simulation of X-ray afterglow error region/galaxy coincidences within 
our 72 GRB fields. Approximately 7 overlaps are expected, while 14 are observed (indicated 
by the arrow). The probability of observing 14 or more overlaps is only 0.5%. 
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Fig. 6. — Observed galaxy properties: (a) galaxy ellipse area versus magnitude, (b) 
ISOAREAF_IMAGE area versus magnitude, (c) ellipticity versus magnitude, (d) galaxy 
ellipse area versus ellipticity and (e) galaxy ISOAREAFJMAGE area versus ellipticity. The 
plot axes limits are chosen to emphasize the afterglow-associated galaxies and exclude some 
of the galaxies in this sample. While not plotted, the galaxies are included in the Monte 
Carlo simulations. The galaxies associated with X-ray afterglows {filled and open squares) 
are a typical sample of all the field galaxies {gray points). The galaxies associated with both 
X-ray and optical afterglows {open squares) differ from the field galaxy population at the 
~ 2a level for all relationships involving magnitude {a,b,c). 



